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We explore the implications of electroweak baryogenesis for future searches for permanent electric 
| dipole moments in the context of the minimal supersymmetric extension of the Standard Model 
(MSSM). From a cosmological standpoint, we point out that regions of parameter space that over- 
■ produce relic lightest supersymmetric particles can be salvaged only by assuming a dilution of the 
& ' particle relic density that makes it compatible with the dark matter density: this dilution must occur 
after dark matter freeze-out, which ordinarily takes place after electroweak baryogenesis, implying the 
same degree of dilution for the generated baryon number density as well. We expand on previous 
studies on the viable MSSM regions for baryogenesis, exploring for the first time an orthogonal slice 
of the relevant parameter space, namely the (tan /?, itia) plane, and the case of non-universal relative 
gaugino-higgsino CP violating phases. The main result of our study is that in all cases lower limits 
on the size of the electric dipole moments exist, and are typically on the same order, or above, the 
expected sensitivity of the next generation of experimental searches, implying that MSSM electroweak 
baryogenesis will be soon conclusively tested. 



1 Introduction 



Explaining the origin of the baryonic matter of the Universe remains an open problem at the interface 
of cosmology with particle and nuclear physics. Consistent values for the ratio Yb of the baryon number 
density {jib) to entropy density (s) have been obtained from measured light element abundances in 
the context of Big Bang Nucleosynthesis (BBN) |T] and acoustic oscillations in the cosmic microwave 
background (CMB) as measured by the Wilkinson Microwave Anisotropy Probe (WMAP) [HE]: 



Assuming that the Yb = at the end of the inflationary epoch, the dynamics of the subsequently 
evolving cosmos would then have generated a non-vanishing baryon asymmetry. As first noted by 
Sakharov, [3J, these dynamics must have included violation of baryon number (B) conservation; viola- 
tion of C and CP conservation; and a departure from equilibrium dynamics assuming CPT invariance. 

Although it is not known when in the history of the universe these ingredients came into play, one 
possibility that can be tested with terrestrial experiments is that baryogenesis occurred during the era 
of electroweak symmetry breaking. In this scenario, which we consider here, a departure from equi- 
librium occurs through a strong, first order electroweak phase transition (EWPT) in which bubbles of 
broken electroweak symmetry form in the symmetric background. Charge asymmetries generated by 
CP-violating interactions at the bubble walls diffuse into the unbroken phase, where B-violating elec- 
troweak sphalerons convert them into non-vanishing baryon number density. The expanding bubbles 
capture this baryon number density by quenching the sphaleron transitions, ultimately leading to a 
relic baryon asymmetry at zero temperature. 

In principle, the Standard Model (SM) of particle physics could have satisfied these "Sakharov 
criteria" as needed for electroweak baryogenesis, particularly as the B-violating sphaleron transitions 
associated with the SU(2) l gauge sector are unsuppressed at temperatures of order 100 GeV and above. 
However, the strength of the SM CP-violating (CPV) interactions are too small to have generated 
sufficiently large charge asymmetries needed to bias the sphalerons into making baryons. In addition, 
the lower bound on the mass of the SM Higgs boson [3] implies that finite temperature SM effective 
scalar potential that governs electroweak symmetry-breaking would not admit a sufficiently strong 
first order EWPT as needed to prevent a "washout" of any baryon asymmetry that might otherwise 
have been produced at high temperatures 0[6]. Consequently, an explanation of the observed Yb 
requires new physics beyond the Standard Model . 

In this paper, we concentrate on the possibility that this new physics involved TeV scale super- 
symmetric interactions in the guise of the minimal supersymmetric Standard Model (MSSM). The 
general motivation for considering TeV scale supersymmetry is well-known, and we refer the reader 
to the vast literature for a discussion [7]. Here, we focus on the minimal scenario, up-dating our 
previous study to take into account several developments (outlined below). The phenomenological 
constraints on MSSM electroweak baryogenesis (EWB) are quite stringent: the existence of a strong 
first order EWPT requires both a relatively light SM-like Higgs scalar and a correspondingly light 
scalar superpartner of the right-handed top quark. Searches for scalars at the Large Hadron Collider 
(LHC) will readily explore the narrow window of available parameter space for these considerations. 
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The strengths of MSSM CPV interactions relevant for EWB are similarly tightly bounded by current 
limits on the permanent electric dipole moments (EDMs) of the electron, neutron, and mercury atom, 
with bounds on the latter having recently improved by a factor of seven [8]. Thus, it is conceiv- 
able that the combination of LHC studies and planned EDM searches having substantially improved 
sensitivities may rule out the minimal version of supersymmetric EWB, leaving other scenarios such 
as non-minimal supersymmetry or leptogenesis as the most viable possibilities. The purpose of the 
present study is to delineate just how close we are to such a situation. 

The new recent developments in supersymmetric EWB include the following results: 

(i) Carena et al. [9] have completed an analysis of the EWPT in the MSSM using a two-loop, finite 
temperature effective potential computed with effective theory to integrate out all but the most 
relevant scalar degrees of freedom. Their results indicate that a sufficiently strong first order 
EWPT can occur if the SM-like Higgs and RH stop masses are less than 125 GeV, with some 
dependence on the cutoff scale of the effective theory. For the EWPT- viable regions of parameter 
space, the electroweak minimum is metastable with a lifetime longer than the age of the universe. 
The presence of the light, RH-stop leads to the existence of a deeper zero-temperature color- 
breaking minimum that can be avoided if the electroweak minimum is deeper at temperatures 
associated with the EWPT. 

(ii) Three of us have recently completed a two-loop computation of the electric dipole moments of 
the electron and neutron in the MSSM, including all of the CPV interactions that are most likely 
to be responsible for EWB in the MSSM [10]. The results indicate that even under the most 
optimistic scenarios, viable MSSM EWB requires that the CPV phases in the gaugino-Higgs- 
Higgsino sector be non- universal, in contrast to the usual universality assumptions. The CPV 
interaction most likely to lead to successful EWB involves the bino and Higgsino fermions 

(iii) It was shown in Ref. |12] that the value of tan j3 could play a decisive role in the transport dy- 
namics that convert the bino-driven Higgsino asymmetry into a left-handed SM fermion number 
asymmety - the quantity that ultimately biases the sphalerons into making baryons. Larger val- 
ues of tan (5 lead to a suppression of Yb and a potentially significant dependence on the masses 
of the RH bottom squarks and tau sleptons (tr), whose Yukawa couplings are enhanced 
with increasing tan/3. For values of this parameter favored by the muon anomalous magnetic 
moment, one expects a suppression of a factor of two or more - and possibly even a sign change 
in Yb - depending on the RH sbottom and stau masses. 

In addition to analyzing the phenomenological implications of some of these developments for 
MSSM EWB and EDM searches, we also take into account another physical consideration not treated 
in previous work, namely the impact of the cold dark matter (CDM) relic density at the time the 
baryons are produced. In some regions of parameter space in which MSSM EWB is naively viable, one 
expects f^cDM due to neutralinos to be larger than the observed relic density. In standard cosmological 
scenarios, these parameter space regions are therefore excluded. They may become viable in non- 
standard scenarios involving late-time reheating that lead to a dilution of the neutralino density due 
to entropy production. This entropy dilution, however, will also lead to a reduction in Yb- We find 
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that the resulting impact on the EWB-allowed MSSM parameter space can be substantial, particularly 
in regions where wino-Higgsino driven EWB might otherwise have been potentially important. The 
impact on the bino-Higgsino driven scenario is less pronounced since in this region of parameter space 
the neutralino relic density typically lies below the observed value so that no entropy dilution rescaling 
occurs. 

Taking all of the aforementioned considerations into account, we find minimum values for \d e \ 
and \d n \ for which MSSM EWB would be viable under the most optimistic scenarios. These values 
lie somewhat below the expected sensitivities of the next generation of EDM searches, though given 
outstanding theoretical uncertainties (discussed below) we cannot make a definitive statement until 
further theoretical progress is achieved. For the interim, we provide this analysis as a phenomenological 
"progress report" that will undoubtedly require further updates in the future. To this end, we illustrate 
the key features of the current status in several figures: (a) a series (Fig. 1-6) that shows the EWB- 
viable regions in the plane of gaugino-Higgsino mass parameters for different values of the CP- violating 
phases and corresponding values of the EDMs. In this series, we show the impact of imposing the 
requirements of entropy dilution needed for consistency with the observed CDM relic density; (b) a 
series of plots (Fig. 7 - 10) in the space of tan/3 and the mass of the CP-odd scalar, uia- 

We discuss in Section [2] the general framework, and we outline the MSSM parameter space com- 
patible with EWB. In this section only we consider a light right handed stop as the driver of a strongly 
first order EW phase transition: in the rest of the paper we drop this assumption, and always consider 
(i) a heavy sfermion sector (including both stops) and (ii) beyond the MSSM physics driving a strongly 
first order EWPT. In Section [3] we discuss the relation between the neutralino relic abundance and 
EWB . For purposes of comparison with our earlier work, we consider here a scenario of universal 
CPV gaugino-Higgsino-Higgs phases without imposing the EDM constraints. The latter are treated 
in Sections H] and [5j Constraints in the tan /3-rriA plane are given in Section |6j We summarize our 
findings and the outstanding theoretical issues in Section [71 

2 MSSM Baryogenesis and EDMs 

For purposes of our analysis, in this section (and only here) we start with the customary assumption 
that the RH stop and SM-like Higgs masses lie in the ranges indicated by the analysis of Ref. [13] , so 
that the universe undergoes a sufficiently strong first order EWPT to prevent baryon number washout 
(we relax this assumption and consider a heavy sfermion sector, including both stops, everywhere 
else in this analysis). Consequently, we consider a scenario under which the leading CPV particle 
asymmetry generation occurs in the gaugino-Higgsino sector. In order to accommodate the lower 
bound on the lightest SM-like Higgs boson, the LH stop mass must be of order one TeV, thereby 
leading to Boltzmann suppression of LH stops and the production of large CPV asymmetries in the 
stop sector. In principle, one could counteract this suppression to some extent through resonant 
CPV stop processes wherein ~ m t i but our assumption of a heavy sfermion sector suppresses 
the importance of this contribution, even if resonant. For relatively large tan/3 leading to enhanced 
bottom and tau Yukawa couplings, one might introduce additional sources associated with these 
scalar quarks, but we leave consideration of this possibility to a future study. The remaining source 
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of CPV asymmetries thus arises in the gaugino-Higgsino sector, for which one may consider relatively 
light superpartner masses (of order a few hundred GeV). The corresponding CPV phases are, after 
appropriate field redfinitions: 

4>j = Arg(/iM ? -6*) , (2) 

where [i is the supersymmetric Higgsino-Higgs mass parmeter; Mj are the SUSY-breaking soft mass 
parameters for the bino (j = 1) and wino (j = 2); and b is a soft SUSY-breaking Higgs mass parameter. 

To obtain Yb as a function of the MSSM parameters, we rely on the work of Ref. |12lll4rfT6] . 
which include the parameter dependence of both the CPV source terms in the quantum Boltzmann 
equations as well as that of the CP-conserving, particle number-changing interactions. Both sets of 
effects are required for a robust prediction: the CPV source terms determine the overall scale of CPV 
Higgsino-Higgs asymmetries that may arise, while the particle number-changing interactions determine 
the efficiency with which these asymmetries convert into the left-handed fermion density that biases 
the electroweak sphalerons. To lowest order in Yukawa and gauge couplings, particle number changing 
processes occur through decays and inverse decays A o B + C; in regions where these reactions are 
kinematically forbidden, A + B <H> C + D scattering reactions dominate. We avoid these parameter 
space regions as there does not yet exist a complete study of the parameter dependence of these 
scattering rates. 

For the CPV sources, we also rely on the computation of Ref. [H], which computed the sources 
to leading non-trivial order in the spacetime varying Higgs background field. This approximation is 
likely to overestimate the source strength 0, so we consider our results for Yb to represent the most 
optimistic scenario, and therefore, the most conservative for statements about the impact of EDMs. 
On-going theoretical work is aimed at resumming the background field while consistently including the 
effects of diffusion ahead of the bubble wall and interactions that would drive the out-of-equilibrium 
plasma toward equilibrium. Once such a consistent framework has been achieved, our analysis will 
likely need to be up-dated. 

We also note that the sources of Ref. [J3] apply to only resonant production of CPV Higgsino 
asymmetries, and we do not consider non-resonant processes that generally require substantially larger 
CPV phases in order to produce the observed Yb (see, e.g. Ref. [T7] )• Resonant EWB occurs when 
the mass parameters of two particles having a trilinear coupling with the Higgs scalar are nearly 
degenerate. In the case of top squarks, for example, resonant asymmetry production occurs when for 
nearly equal, and light enough LH and RH stop masses, a situation precluded by the requirements of 
a strong first order EWPT and the bound on the lightest Higgs mass as noted above. For the gauginos 
and Higgsinos, resonant EWB arises when /i ~ M\^, - the scenario on which we concentrate herqj. 

Within this context, the overall source strength depends crucially on both the CPV phases and 
A/3, the change in tan" 1 v u /vd from the bubble wall exterior to its interior with v u and Vd being the 
spacetime-dependent vacuum expectation values (or background fields) associated with the neutral up- 
and down-type Higgs scalars, respectively. The authors of Ref. [20] showed that A/3 decreases with 
increasing m^, and we take this m^-dependence into account in our analysis of the parameter space. 
Doing so is not entirely consistent with the phase transition analysis of Ref. [13], which considered 

^ee, e.g. Refs [IMS] 

2 In our conventions, the Higgs and Higgsino fields have been rotated so that fj, is real and positive 
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m,A larger than ~ 10 TeV. We expect that a two-loop EWPT analysis allowing for a dynamically 
active, light CP-odd Higgs would not change the allowed ranges of and m/j substantially, and we 
therefore consider values of tua below one TeV as needed for unsuppressed, resonant, CPV sources. 

The dependence of Yb on tan/3 remains a topic of on-going research. To our knowledge, there 
exist no results in the literature, for the tan ^-dependence of A/3. On the other hand, a previously 
overlooked tan /3-dependence of transport dynamics has been observed in Refs. [T2]. Earlier work had 
assumed that the CP-conserving, particle number changing processes were dominated by stop quark 
(squark) Yukawa interactions, given the large value of the top Yukawa coupling relative to that of 
the other (s)fermions. For moderate tan /3, however, the (s)bottom and (s)tau Yukawa couplings are 
relatively enhanced, and their impact on the conversion of a Higgs- Higgsino asymmetry into the net 
left-handed fermion density {n\ c ^) that couples to electroweak sphalerons can be substantial. In what 
follows, we include this tan /3- dependence by re-scaling Yb according to Fig. 1 of Ref. |16j . that shows 
the ratio of Yb computed including the bottom and tau Yukawa couplings to the baryon asymmetry 
computed in their absence as a function of tan /3. For purposes of this analysis, we work in the limit 
of heavy tr , but note that - as shown in Fig. 1 of Ref. [16] - the presence of a light tr can lead to 
a suppression of Yb- We also work in the region of parameter space wherein >> to avoid 
further suppression of Yb or even a sign change that can occur for mg^ 1z m i R ■ 

For the computation of bounds from EDMs, we draw on the two-loop computations of Ref. |1(J| . 
Limits from one-loop EDMs generally imply values of | sin (ftj | that are too small to accommodate the 
observed value of Yb, even in the regions of resonant EWB. To circumvent these bounds, we work in 
the limit of heavy first generation sfermions, thereby suppressing one-loop EDMs without quenching 
Yb via Boltzmann suppression of stops, gauginos, and Higgsinos. In this limit, both the elementary 
fermion EDMs and the chromo-EDMs of quarks are suppressed. As a result, constraints from the 
i99jjg a ^ om j c EDM limit are minimal, as the 199 Hg atomic EDM is dominated by the chromo-EDMs 
in the MSSM. In contrast, significant bounds from the two-loop EDMs of the electron and quarks 
(contributing to d n ) can occur, with the most stringent arising from graphs involving closed chargino 
and neutralino loops and the exchange of two bosons as in Fig. 1 of Ref. [ID]. To analyze these 
constraints, we rely on the work of Ref. [ID] and consider two cases: (a) universal gaugino phases, 
wherein <p\ = 02, and (b) non- universal phases, for which 0i and 02 are allowed to differ. For the 
second scenario, the sensitivity of the electron and neutron EDMs to 02 is roughly fifty times stronger 
than the sensitivity to 0i. 

In summary, the most optimistic scenario leaving the largest window for viable MSSM EWB 
occurs in the following region of parameter space: a light right-handed stop with m^ R ^5 125 GeV; 
the remaining sfermion masses being heavy, rrij ^ 1 TeV; light gauginos and Higgsinos with \i ~ M\\ 
non-universal phases: 0i ^ 02 j a relatively light pseudoscalar Higgs and moderate to small tan/3. In 
what follows, we illustrate the impact of relaxing these assumptions. 

The MSSM parameter space compatible with EWB is severely constrained by the requirement 
that the lightest neutralino be the lightest supersymmetric particle (LSP), since a light RH stop with 
nii R < 125 GeV forces the LSP mass to also be below that valuqf|. However, the requirement of a 
light RH stop can be relaxed in several minimal extensions to the MSSM that include for instance 

3 A caveat is the assumption of a gravitino LSP, but we do not consider this possibility here. 
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an extended Higgs sector. One large class of such models are those that include a gauge singlet 
with couplings to the SU(2) Higgs doublets [21J. Other viable beyond-the-MSSM setups featuring a 
strongly first order EWPT also include models with a fourth fermion generation [22] and top-flavor 
models [23]; models with higher dimensional scalar representations of SU(2) (which might lead to a 
multistep EWPT [23]) provide yet another promising such framework [25]. In several of these models, 
the physical mechanism that drives a strongly first order EWPT is decoupled from EWB as the source 
of the BAU. As such, in what follows we will assume that the RH stop be at the TeV scale, and 
that a physical process that does not interfere with EWB drives a strongly first order EWPT. This 
assumption will allow us to explore wider regions of the parameter space relevant to EWB, although 
we shall also show the regions that would correspond to an LSP lighter than 125 GeV (dashed blue 
lines in fig. Q]and[5]). From now on, in this manuscript we take all sfermions, including the RH stop, 
to be heavy, m? ^ 1 TeV. 

3 Baryogenesis and Neutralino Relic (Over-) Abundance 

If the entropy density of the Universe is not a constant after the freezeout of the net baryon density, 
the resulting baryon asymmetry today will in general be affected. While episodes such as very low- 
scale inflation (as low as the electroweak scale) have been envisioned |26H28j . this is not the standard 
paradigm. However, an entropy injection episode might actually be needed in a successful theory of 
supersymmetric electroweak baryogenesis. We provide here an example and study the implications for 
the underlying supersymmetric theory. 

In the framework we consider here, where i?-parity is conserved, the LSP is stable. Given null 
results in searches for a stable massive charged or strongly-interacting particle @, the LSP must be 
neutral, and its density today must be compatible with (namely, at least not exceed) the universal 
matter density to avoid over-closing the Universe. As such, a relic, neutral LSP x must have abundance 
Q x h (in units of the critical density) less or equal to the dark matter density Odm- If this is not the 
case, either the underlying supersymmetric theory is ruled out or a mechanism must exist that dilutes 
away the excessive x relic abundance. 

Weakly interacting massive particles, such as the lightest neutralino in MSSM scenarios where 
it is the LSP, undergo a rather universal (i.e. independent of the details of masses and interac- 
tions) freezeout process from the thermal bath in the early universe that occurs around temperatures 
^freeze-out ~ m x/(^0 ^5). In models of successful electroweak baryogenesis, where the LSP cannot 
be heavier than ~ 1 TeV [30], the freezeout temperature is thus always below 50 GeV. Therefore, 
freeze-out occurs typically at much lower temperatures than those relevant for the electroweak phase 
transition. In turn, this implies that if an entropy injection episode is responsible for the dilution of 
an excessive thermal relic density of neutralinos, the same entropy injection episode will also dilute 
away the baryon number density produced at the electroweak phase transition. 

We are therefore left with two possibilities: either models with an excessive neutralino thermal 
relic abundance are ruled out, or they are salvaged by a modified cosmological expansion that inflates 

4 In particular, the number density of a metastable negatively charged electroweak-scale particle \~ relative to entropy 
is constrained to be smaller than n x -/s < 3 x 10~ 17 , since it leads to catalytic enhancement of the 6 Li production. [29) 
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them away □• However, if this second possibility is realized, the net baryon asymmetry will also be 
diluted away, requiring a larger asymmetry to be produced at the electroweak phase transition. It is 
this latter possibility that we explore in this section. 



5 Interestingly, such a mechanism can be built in the electroweak phase transition itself, see e.g. the recent analysis of 
Ref. EH 
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Figure 1: Iso-level curves for the lightest neutralino relic abundance (upper panels) and for the gaugino-higgsino 
CP V phase producing the central value of the BAU, without (central panels) and with (lower panels) entropy 
rescaling for models with over-abundant relic neutralinos. All panels refer to the (M\, fi) plane, with itla = 150 
Ge V in the left panels, and tua = 500 Ge V for those to the right. The red shaded region is excluded by the non- 
observation of light neutralino pairs at LEP. The green bands correspond to a neutralino relic density consistent 
with the WMAP results. The blue dashed lines in the upper panels indicate a lightest neutralino mass of 125 
GeV. 
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To illustrate the effect of enforcing a cosmologically acceptable density of relic neutralinos, we 
consider in the upper two panels of Fig. [T] the thermal relic density of neutralinos in the (Mi, fx) plane, 
for a scenario with heavy sfermions, tan/3 = 10, M2 = 2 x Mi = M3/3 and txia = 150 GeV (top left 
panel) and 500 GeV (top right panel). The green line corresponds to a thermal relic abundance in 
accord with the inferred density of cosmological dark matter. In the lower right corners of these two 
plots, neutralinos only constitute a fraction of the cosmological dark matter (unless some non-thermal 
production mechanism is operative). Although this would mean one needs to postulate an additional 
particle to provide the rest of the dark matter density, an under-abundance of thermal neutralinos is 
cosmologically viable, and there is no need to rescale the baryon asymmetry that gets produced in 
those regions of the parameter space. The blue dashed lines in the upper panels indicate a lightest 
neutralino mass of 125 GeV. 

In the lower right corner fj, < Mi, and the lightest neutralino is higgsino-like, and it efficiently 
pair-annihilates into W + W~ ad ZZ pairs. Further contribution to the efficient pair- annihilation 
processes in the early Universe also stems from the co-annihilation of the quasi-degenerate next-to- 
lightest higgsino-like neutralino and chargino. In the upper right portion of the plots, the increasing 
bino-like content of the lightest neutralino mass eigenstate suppresses annihilation to the SU(2) gauge 
bosons, and only if a resonant annihilation channel occurs (e.g. for m x ~ myi/2 in the left plot) can 
the relic abundance be compatible with the upper limit to the cosmological dark matter density. In 
that region, therefore, we will need to rescale the neutralino density by at least a factor O^/Odmj the 
former factor being the thermal neutralino relic density in units of the critical density. 

Without accounting for such rescaling, the middle two panels in Fig. 1 show the values of sin0 
one would need to produce the observed baryon asymmetry of the Universe. Each contour of the 
double funnel-like region corresponds to a given value of shi(/>, where we take (p\ = $2 = 4> under 
the assumption of phase universality. Going from the outermost funnel to the innermost, successive 
contours correspond to the values of s'm<j) indicated in the panels "isolevels" ) . The behavior of the 
contours is consistent with the results given in Fig. 1 of our earlier work |30j. Increasing the value of 
rriA suppresses the produced net baryon asymmetry, and therefore the viable regions (those inside the 
funnels corresponding to the resonant regions where (i ~ Mi or ~ M2) shrink. 

Accounting for the neutralino number density dilution and transferring the dilution to the baryon 
number asymmetry as well leads to a requirement of larger sm<p (i.e. a larger initial net baryon number 
density) in regions where Q x > ^dm- This feature is illustrated in the bottom two panels of Fig. 1. 
There, we find a suppression of the M2 ~ fi funnel, where the lightest neutralino is dominantly bino- 
like. For large values of txia (lower left panel), said funnel essentially disappears, leaving bino-driven 
electroweak baryogenesis [TT] as the only available option. 

In the following section we further explore the phenomenological consequences of entertaining a 
viable supersymmetric model from the standpoint of cosmology, and we derive lower limits on the 
electron and neutron EDM in the present scenario. 
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Figure 2: Curves at constant CPV gaugino-higgsino phase, on the (Mi, fx) plane, at vtla = 300 GeV and 
M2 = 2M\, with (right) and without (left) entropy rescaling for overabundance relic neutralinos models. The 
grey region does not produce a large enough BA U through resonant processes, the red region is excluded by LEP 
searches for the lightest chargino, and the green lines indicate a neutralino relic abundance equal to the cold dark 
matter abundance in a standard cosmological setup. 

4 A Lower Limit on the Electron and Neutron EDM 

Fig. [21 left, reiterates what is shown in the Fig. [H but this time for a heavy Higgs sector set to 
rriA = 300 GeV. As in the previous plots, we again assume here M2 = 2M±. The different shape in the 
contours depend on the assumed values for m\. Again, the green lines indicate a thermal neutralino 
relic density in accord with the inferred cosmological dark matter density. The funnel where the 
lightest neutralino resonantly annihilates through the s-channel on-shell exchange of one of the heavy 
neutral Higgses appears near M\ ~ m^/2 = 150 GeV. The left panel shows contours of constant values 
for the gaugino-higgsino phase needed to have the right amount of baryon asymmetry, neglecting the 
overabundance of neutralinos in the upper- left part of the parameter space. As indicated by the right 
panel, the effect of the required entropy dilution shrinks the available regions for those values of (Mi, [i) 
where f2 x > SIdm- 

Focusing now on the physically motivated case of a re-scaled net baryon asymmetry, Fig. [3] shows 
lines of constant values for the electron (left) and for the neutron (right) EDMs. As one would expect, 
the largest values occur where the CP violating phases are largest, i.e. near the outmost boundaries 
of the funnel regions, while smaller values are obtained in the central part of the funnels: there, 
the resonant source term allows for a lower CP violating phase, resulting in smaller electric dipole 
moments. The only viable regions are therefore those inside the blue contours. In summary, we 
find for the slice of parameter space we consider that there is a lower limit to the electric dipole 
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Figure 3: Curves of constant values for the electron (left) and for the neutron (right) electric dipole moment, 
on the same plane as Fig. [H right. The regions outside the blue contours are excluded by current experimental 
limits. 

moment of the electron around 0.4 x 10~ 27 e cm, roughly a factor 4 below the current experimental 
limit; the analogous limit for the neutron EDM sits at 10 -26 e cm, i.e. only a factor 3 below current 
sensitivity. We emphasize that these lower limits pertain only to scenarios involving universality of 
phases: Arg(/iMi&*) = Arg(/iM2&*) = (j)^. In section 5, we will relax this unnecessary assumption. 

Up to this point, we have considered only the (complete set [TO]) two- loop contributions to the 
EDMs, since one-loop contributions are negligible for sufficiently heavy 1st generation sfermion masses. 
The following Fig. [4] addresses the question of how small the sfermion mass scale must be for one loop 
contributions to become dominant over the two- loop ones. The left and right panels of Fig. 0] focus on 
the case of the EDM of the electron and of the neutron, respectively. We set all CP violating phases 
to zero with the sole exception of the relative higgsino-gaugino phase, as in the rest of this study. The 
sfermion scale we refer to here alludes naturally to the relevant sfermions, i.e. first generation sleptons 
and squarks. Roughly, the size of the CP violating phase factors out from both the one- and two- loop 
contributions. The larger the values of (Mi,//), the more suppressed the two- loop contributions: as 
a result, the sfermion mass scale that would produce a comparable EDM can also be larger, as the 
iso-level contour labels show. Fixing the sfermion mass scale to e.g. 6 TeV, points to the upper-right 
of the lines have a larger one-loop contribution than the two-loop one, while points to the lower-left 
side of the lines have a larger two-loop contribution. 

In summary, we find that the one-loop contributions to the electron EDM start to be significant 
when the mass scale is around 4-7 TeV, while squarks as heavy as 10 TeV can give significant contri- 
butions at the one-loop level to the neutron EDM. This implies, in particular, that the neutron EDM 
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Figure 4: Curves indicating the values of a common sfermion mass scale such that one and two loop contri- 
butions to the electron (left) and to the neutron (right) EDM are equal, on the same plane as Fig. [H right. 
At a given sfermion mass scale indicated on the lines, points inside the funnel region lying closer to the origin 
correspond to two-loop contributions that are larger in magnitude than one-loop contributions. The converse 
holds for points lying closer to the funnel tips. 

is overall more sensitive to one-loop contributions than the electron EDM 

5 EDM in the bino-driven electroweak baryogenesis scenario 

As pointed out by some of us in Ref. [11], it is possible for MSSM baryogenesis to remain viable even 
if the next generation of EDM searches yield null results. Theoretically, this scenario requires that 
one relaxes the assumption of universality between the relative bino-higgsino and wino-higgsino phase, 
i.e. 4>i 7^ 4>2 and work in the limit of heavy 1st generation sfermion masses. In this way, as long as 
<f>2 is sufficiently small, EDMs are suppressed, while bino-driven electroweak baryogenesis, efficient if 
the resonant condition /i ~ Mi is approximately valid, can proceed with large bino-higgsino phases 
4>i. The latter, as shown and explained in detail in Ref. |11| . does not lead to excessively large 
contributions to the EDMs, even if | sin | is of order one. 

In the bino-driven portion of the parameter space, Sl x < Hdm, and no entropy dilution is needed 
to wash out excess relic neutralinos. We illustrate this point in the left panel of Fig. 5, where we 
assume ^2 = and show the values of 4>\ for which we generate the observed baryon asymmetry. Since 
resonant bino-driven EWB is most efficient in the region for fi ~ Mi, it is useful to re-cast our results 
in terms of Mi and the relative mass-splitting (// — M\)/M\. In the right panel of Fig. 5, we show the 
bino-driven region (i.e. the region where [i ~ Mi) in the (Mi, (jj, — Mi)/M\) plane. As the contours 
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Figure 5: Curves indicating constant values of the relative bino-higgsino phase <f>i that produces the right BAU 
through resonant processes, on the same plane as fig. [H (left), and on the plane defined by M\ and by the relative 
bino-higgsino mass splitting (p — M\)/Mi. The grey region does not produce a large enough BAU through 
resonant processes, the red region is excluded by LEP searches for the lightest chargino. In the left panel, 
the green lines indicate a neutralino relic abundance equal to the cold dark matter abundance in a standard 
cosmological setup. In the right panel, the blue dashed line indicates a lightest neutralino mass of 125 GeV. 

for successively smaller values of sin0 indicate, a lighter overall mass scale and smaller relative mass 
splitting are needed to reproduce the observed Yb as sin^ is decreased. 

We employ the same parameter space (the (Mi, (/i — M\)/M{) plane) to illustrate that even in 
the case of bino-driven EWB we obtain an absolute lower limit on the size of the electron and the 
neutron EDM that would be consistent with this scenario. Here we give contours of constant d e (left 
panel) and d n (right panel), rather than constant sin^i, consistent with Yb- The contours in the left 
panel of Fig. O left, indicate that the minimal electron EDM is around 10 -29 e cm, roughly a factor 
200 below present limits. The projected sensitivity improvement to cover this slice of the bino-driven 
electroweak baryogenesis scenario with neutron EDM searches is instead of around a factor 100, as 
shown in the right panel. 

6 The (tan /3, niA) Plane 

Thus far, our exploration of the supersymmetric parameter space compatible with EWB has focused 
largely on the parameter space defined by the masses of those particles directly entering the source 
terms in the production of the baryon asymmetry at the electroweak phase transition. We now turn our 
attention to two other parameters that play a crucial role in the dynamics of baryon number generation 
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Figure 6: Curves indicating constant values for the electron (left) and for the neutron (right) electric dipole 
moment, on the same plane as fig. [5[ right. The CPV phase <pi is set to the value giving the right BAU, while 
4>2 (the relative wino-higgsino phase) is set to 0. 



at electroweak temperatures: the ratio of the vacuum expectation values of the neutral componets of 
the two SU(2) Higgs doublets, tan/?, and the mass of the CP-odd Higgs, m,A- The former enters in 
the calculation of Yg via not only directly in the sources via the various particle couplings, including 
Yukawa couplings, but also via the effects described in Ref. [T2"l[l"6] . In the latter work, it was observed 
that for moderate values of tan/3, the timescale for b (squark) and tau (s)lepton Yukawa interactions 
can be shorter than the timescale associated with diffusion ahead at the bubble wall. Consequently, for 
sufficiently light right-handed sbottoms and staus, the dynamics of these sfermions can substantially 
alter the transport dynamics that govern the conversion of Higgs-Higgsino CP-violating asymmetries 
into the (s) quark sector asymmetries. The mass parameter rriA also enters in the calculation of Yg 
through the parameter A/3, to which Yg is linearly proportional. Here, we treat the dependence of Yg 
on iua according to the two- loop results presented in [2U] , and the dependence on tan /3 according to 
the full treatment outlined in [12] [§. 

So far, the (tan [3, ttia) parameter space for this problem has not been extensively explored. For 
illustrative purposes, we adopt here the universal phase assumption and we show that there can be 
strong bounds in this space that have not been emphasized in previous work. At the end of this 
analysis, we discuss briefly our expectations for the case of non-universal phases. 

In Fig. 7, we show the constraints on the (tan/3, itia) parameter space from various phenomenolog- 
ical and cosmological considerations. We set the gaugino and higgsino mass parameters to M± = 145 

6 Note that we do not consider the effects of varying ttia on the strength of the electroweak phase transition since in 
the latest work of Ref. [9] , the pseudoscalar A was integrated out of the effective theory. 
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m =145 GeV, m =290 GeV, li=300 GeV 




Figure 7: The (tan/3, tua) plane: the orange region is ruled out by LEP limits on the lightest Higgs, while the 
grey region is ruled out by excessive supersymmetric contributions to the inclusive b — > sj branching ratio and 
the yellow region is ruled out by bounds on the BR(B U — > tv t ). In the blue region at the top the generated BAU 
is too small, while the green contours indicate the 2-a preferred range for relic neutralinos to be the dark matter. 
We also show contours of constant values for the neutralino relic abundance (green dashed lines) and for the 
value of the gaugino-higgsino CPV phase sin^>„ (solid blue lines) such that electroweak baryogenesis successfully 
accounts for the BA U. 



GeV, M2 = 290 GeV and \x = 300 GeV. As discussed above, we also assume a common gaugino- 
higgsino phase <j) = Arg(/uMj6*), j = 1,2. In the orange region to the left and bottom of the plane, 
the Higgs mass is lower than the limit set by LEP-II of 114.4 GeV (this limit applies for the heavy 
sfermion sector we consider in the present analysis). In addition, in the grey area the supersymmetric 
contribution to the b — > 57 decay exceeds the experimentally allowed value. 

At large values of tan/3, a significant constraint stems from MSSM contributions to the B u — > tv t 
decay mode [3"2"ll3"3"] . This branching ratio, in turn, depends on SUSY-QCD corrections to the charged 
Higgs boson coupling to fermions. These corrections are asymptotically small in the limit of heavy 
sfermions, and, for our choices of parameters, the relevant quantity eo in Eq. (6) of Ref. |32j is always 
smaller than 0.01. Given that in the present analysis we mostly focus on the limit of heavy sfermions 
to avoid one loop contributions to the EDMs, we set eo = 0, which entails slightly more stringent 
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constraints than for a small, finite positive value. Imposing the 2-a limit to the ratio of the MSSM 
to SM contribution 0.53 < R TUt < 2.03 [22], two disconnected regions of the (tan/3, m^) parameter 
space are ruled out. As shown in Ref. [33], the portion of parameter space between these disconnected 
regions is actually also ruled out by constraints on MSSM contributions to the ratio of K — > [iv^ and 
7r — > iiv^ branching ratios. For the range of tan /3 shown in Figs. 7-10, the latter region also lies within 
the b — > 57 exclusion region. We shade in yellow the remaining excluded parameter space implied by 
the B u — > tu t constraints. 

The green dashed lines correspond to various, constant values of the lightest neutralino relic abun- 
dance Q x h 2 (as indicated by the labels). Within the thick green bands, the thermal neutralino relic 
abundance is within the 2a WMAP range for the average cold dark matter density. Finally, the blue 
lines indicate which value of the phase <f> is needed for the particular point in parameter space to pro- 
duce the observed Yb- For these lines, we do not include the dilution of over-abundance neutralinos. 
In the upper region of the plot, shaded in cyan, no resonant electroweak baryogenesis can produce 
enough baryon asymmetry, even neglecting the effect of entropy dilution: this therefore sets, for the 
particular values of gaugino and higgsino parameters we picked here, an upper limit of mA < 800 
GeV. 

In Fig. as in the remainder of this paper unless otherwise specified, we have assumed heavy 
squarks and we have neglected the effect of down- type sfermions (effectively assuming heavy sbottoms 
and staus). We study in fig. [8] the impact of two effects: (i) rescaling of the net baryon asymmetry 
for parameter space points with a relic thermal neutralino abundance in excess of the dark matter 
density, and (ii) factoring in finite masses for third generation sfermions. The lines in Fig. [8] represent 
iso- level contours of (corresponding to the observed Yb) under different assumptions. The red lines 
correspond to the limit of super-heavy down-type sfermions, after rescaling for over-dense neutralino 
relic densities. As evident, the parameter space shrinks with respect to the case without density 
rescalling (Fig [7]), and for reasonably low values of tan (3, the maximal value ttia can take is below 600 
GeV. Note that in contrast to the isolevel contours of Fig. [7J those of Fig. [5] increase monotonically 
with tan /3. 

The black lines give the isolevel curves when the effect of bottom and top Yukawa interactions 
are included in the transport dynamics, for right-handed sbottom and stau masses equal to 1 TeV. 
Inclusion of these RH sfermions suppresses the baryon asymmetry, implying a need for a lighter m\ 
for a given value of sin (p^ . The impact of the RH third generation sfermions is even more pronounced 
when they are relatively light. To illustrate, we show in green the isolevel contours for 300 GeV 
6r and tr masses. The appearance of "elbows" in these contours arises because the suppression of 
Yb in the presence of light 6r and tr grows with tan/3. For the relatively flat portion of the green 
curves at lower tan /3, the suppression due to the bu and tr essentially compensates for the monotonic 
increase that would otherwise occur when these sfermions are heavy (the black curves) . Eventually, the 
light RH sfermion suppression takes over, pushing the value of Yb below the experimental value and 
necessitating a smaller value of txia as needed to enhance the CPV source in the transport equations. 
This transition point corresponds to the elbow in a given isolevel curve. Interestingly, the resulting 
parameter space allowed by the requirement of successful EWB for 300 GeV sfermions is now bounded, 
with tan /3 < 20 and tua < 350 GeV. Note that we have not yet imposed EDM constraints. 
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Figure 8: The same plane as in fig. but taking into account the effect of relic neutralino dilution in models 
over-producing dark matter on the CPV phases needed to achieve successful electroweak baryogenesis. The color 
code is the same as in fig. [^] We indicate with different colors lines corresponding to different assumptions 
on the mass of the Sbottom and of the Stau sfermions. Namely, our default assumption of heavy sfermions is 
indicated in black, while light (300 GeV) sfermions correspond to the green line. Neglecting the effect of the 
Yukawa interactions involving sbottoms and staus produces the lines indicated in red. 

We now investigate the size of the permanent EDMs of the electron (Fig [9|) and of the neutron 
(Fig [TOl) over the same parameter space, taking into account the relic density rescaling, but still 
assuming super-heavy third generation sfermions. Our results are thus conservative here, given that 
including the effect of lighter sfermions would suppress Yb and force larger CP violating phases. 

The shape of the contours in the figures illustrates how the size of the two-loop EDMs depends on 
both the size of uia (larger values suppress the two loop contributions where the heavy Higgs sector 
enters) and on tan [3. In addition, the EDMs also linearly depend on the size of the CP violating 
phases (that partly compensate for the suppression induced by larger values of m^)- As a result, 
the parameter space is bounded, and there is an absolute lower limit to the size of the EDMs on 
this slice of parameter space as well. Specifically, for the choice of parameters we make here, the 
electron EDM exceeds the current experimental limit for itia > 330 GeV and tan f3 > 15, with even 
more stringent constraints from the neutron EDM. The smallest the electron EDM can get on this 
slice of parameter space is a few xlO -28 e cm, less than one order of magnitude below the current 
experimental sensitivity. The neutron EDM (fig. [T0"|) can at most be a factor of a few below the current 
sensitivity. 
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Figure 9: Contours of constant values for the two-loop electron EDM, on the same plane as the one shown 
in fig. \$ (i.e including entropy rescaling). The blue line indicates the current experimental limit: the viable 
parameter space lies below that curve. 

In short, the study of the (tan j3, tua) parameter space enforces the conclusion we obtained from 
the study of the orthogonal gaugino-higgsino parameter space: the size of the electron and of the 
neutron EDM, even only at the two-loop level, is bounded from below, and EWB will be thoroughly 
tested in the next generation of EDM search experiments. For the case of non- universal phases, we 
expect the EDM constraints on the (tan /3, itia) plane to be weaker, given the reduced sensitivity of 
EDMs to the bino phase (assuming it is the one responsible for resonant EWB). With an improvement 
of 0(100) in EDM sensitivity, however, we expect the restrictions on the parameter space in this plane 
will be similar to those shown for the universal phase case. 

7 Discussion and Conclusions 

In this paper we continued our program of investigating the phenomenology of supersymmetric models 
with successful electroweak baryogenesis. Our primary findings are: 

1. Regions of parameter space with over abundant thermal relic neutralino production are either 
cosmologically ruled out or must occur in conjunction with a late-time reheating episode. In the 
latter instance, the consequent dilution of both the neutralino relic density and baryon number 
asymmetry implies that the baryon number density at the EWPT must have been larger. These 
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Figure 10: As in fig. [P| this time for the two-loop neutron electric dipole moment. 

considerations lead to more demanding requirements on the MSSM parameters in the domain 
of resonant wino-Higgsino driven baryogenesis. As a result, we find that the parameter space 
for resonant wino-Higgsino driven baryogenesis is significantly more constrained than previously 
realized. In contrast, the requirements for successful bino-Higgsino driven baryogenesis are 
largely unaffected. 

2. We explored for the first time the parameter space defined by the pair tua and tan f3. We find 
that a general upper limit exists to the maximal value the heavy Higgs mass scale can take, 
and this limit is even stronger if relic density dilution effects are taken into account. If third 
generation sfermions have a low mass scale, then an upper limit to tan (3 also emerges due role 
played by the third generation (s)fermions in the transport dynamics. 

3. General lower limits exist on both the electron and the neutron EDMs, which are even stronger 
when the dilution effects are taken into account. Assuming phase universality, these lower 
bounds lie below the current experimental bounds by less than one order of magnitude. The 
non-observation of electron and neutron EDMs in experiments having ten times better sensitivity 
would likely rule out MSSM baryogenesis under the assumption of phase universality. These 
prospects can be seen from both the accessible regions of the (fi,Mi^) parameter space and in 
the (tan /3,Ma) plane. 

4. If one relaxes the assumption of phase universality, then the expected lower bounds on \d e \ and 



19 



\d n \ consistent with the scenario of resonant, bino-driven MSSM baryogenesis are roughly ten 
times smaller than we find under the assumption of phase universality. Because the strength 
of the CP-violating sources used in our analysis is likely to over-predict the strength of these 
sources that would likely result from a fully consistent, all-orders Higgs vev resummation, we 
expect that the resonant bino-driven scenario will be testable with the next generation electron 
and neutron EDM searches that aim for two orders of magnitude improved sensitivity. 

The general picture that emerges from our present phenomenological update is one where MSSM 
electroweak baryogenesis is a predictive framework that will be conclusively tested in the near future. 
Indeed, the parameter space we combed in this analysis is devoid of regions where the predictions for 
specifically EDMs are unobservably small. Thus, the planned improvements in experimental sensitivity 
will put the entire setup under thorough scrutiny. Although there remain important open theoretical 
questions pertaining to the transport dynamics, such as the tan /3-dependence of the bubble wall 
profiles and consequences of a consistent resummation of the background field for the CP-violating 
sources, we expect that our general conclusions regarding the testability of MSSM baryogenesis will 
continue to hold after these issues are conclusively addressed. 
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